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Abstract 

The general theory of neutron scattering is presented, valid for the whole domain of slow 
neutrons from thermal to ultracold. Particular attention is given to multiple scattering which 
is the dominant process for ultracold neutrons (UCN). For thermal and cold neutrons, when 
the multiple scattering in the target can be neglected, the cross section is reduced to the known 
value. A new expression for inelastic scattering cross section for UCN is proposed. Dynamical 
processes in the target are taken into account and their influence on inelastic scattering of UCN 
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is analyzed. 
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1 Introduction 

Thermal and cold neutrons with wave length 0.03 nm < A < 1 nm is an important tool for 
investigation of condensed matter. Due to the absence of charge and considerably weak inter- 
action with electrons and nuclei incident neutron wave goes deep into target almost without 
distortion and coherently influences on all atoms of the target. The whole specific features 
of the matter (crystalline and magnetic structure etc.) show themselves in interference of the 
secondary scattered waves. 

Theory of neutron-substance interaction for thermal and cold neutrons is well established 
(see, e.g., ^, ||]). It is based on the use of Fermi pseudopotential 

V(r) = ^V v (r - R„) = Yl — M(r - R v ). (1) 

Here r and are the position vectors for neutron with mass m and v-th nucleus with mass M v , 
m u = mM v /{m + M u ) is their reduced mass, a u is the amplitude of neutron scattering on free 
nucleus, connected with the scattering length a v and impact momentum in the center-of-mass 
system k u , in linear on k v a v approximation (valid for slow neutrons), by 

a u = a u {l - ik u a u ). (2) 

For thermal and cold neutrons rescattering of secondary waves is unimportant and one may 
use Born approximation that gives for double differential cross section per one target nucleus 

d 2 a k' 

— - — : K b v'X(w', k, uj). (3) 



dQduj 2nNk 
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Here x{ vv \ K > w ) is the Fourier transform 

+00 

x(vv', K ,lj)= f x(w',K,t)e iuJt dt (4) 



of correlation function 

X {vi/,K,t) = (i\e- iKA »®e iKA »'W\i), (5) 

k, = k — k' and to = e — e' are the neutron momentum and energy transfers. The quantity 
b v = (m/m^di, is called scattering amplitude on bound nucleus, and R,(i) is time dependent 
Heisenberg operator of nuclear position. 

It is easy to show that rescattering of the secondary waves may be indeed neglected for 
thermal and cold neutrons. Coherent summation of the secondary waves may be formed only 
from the volume ~ A 3 . The sum of all amplitudes from this domain is by the order of magnitude 
~ nA 3 (fe/A), where n is the density of the nuclei (~ 10 22 cm" 3 ), and b ~ 1CP 12 cm. So far as 

n&A 2 < 1, i.e. A < 100 nm, (6) 

rescattering can be neglected. For ultracold neutrons (UCN), when A > 10 nm, rescattering 
of neutron wave in media is very essential, and for the neutron wave vector k which satisfies 
k 2 < Annb, rescattering becomes the dominant process and results in the total reflection from 
the surface of the target (of cause, for positive b). Thus, Born approximation in general, and 
the cross section @ in particular, can not be used for UCN. 

If the Born approximation is unjustified then one should start from an exact Schrodinger 
equation for the scattering problem. As the first step one may use a target model with fixed 
(immovable) nuclei and consider integral equation 



* k (r) = e lkr - —2 / , -V(r')V h (r')d 3 r'. (7) 



With a formal use of the Fermi pseudo-potential (|T|) equation (]?]) transforms into 

ifc|r-R„| 

*k(r) = e ikr - J2 K jz t~> i *k(R„)- (8) 



V 



R 



The quantity ^(Rz/) seems to have the meaning of the neutron wave amplitude on the i^-th 
nucleus. This wave is composed now from the incident wave and all reflected waves, so its value 
is not known in advance and consistent equations have to be formulated for these quantities. 
It is impossible to get the mentioned equation from (||) just by substitution there for neutron 
position r = R^, due to infinity in diagonal term of the right-hand side. So, it was in fact 
postulated that a proper equation may be obtained just by throwing away the diagonal term 
(self-scattering) : 

ifc|R„-R„,| 

* k (R„) = ^ ~ E 6 -'T5 s-r*k(Rv). (9) 

v , +v IR^-Rz/I 

The equations (||) and (||) are today the basis for the whole theory of ultracold neutron 
interaction with matter (see, e.g., [|], ||, 0]). From @ one can get an effective repulsive 
(optical) potential on the condensed matter surface, so the neutron wave with the energy below 
the threshold is exponentially decreasing deep into target. 

In general the problem of neutron interaction with matter can be analyzed in the framework 
of Multiple Scattering Theory (MST) §, % 0, MST deals with interaction of a projectile 



with many-body target. In this theory a formal solution of many-body problem takes the form 

# = f°ir 1 VtX (10) 
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with functions ^ v defined by the set of linear equations 

^ = *° + D- 1 £ iu'Vu>, (11) 

where operators t v (i-matrixes) are linked with potentials V v by 

i v = V V + VuD'%. (12) 

Here is the wave function for noninteracting projectile and target, and is the "Green 
function" (see details in the next section). 

There is a wide spread opinion in literature that (||) and @ are just the equations of 
MST (|io| ) and (11) for fixed nuclei target. In fact, they have similar structure but with some 
modeled t v . Consistent derivation of (10) and (|TT| ) analogies for fixed nuclei target with realistic 



potentials was done in [12]. In any way, target model with fixed nuclei can be applied only to 
elastic scattering. While UCN escape from vessels, that is studied for many years, is due mainly 
to inelastic scattering. Now quasielastic scattering of UCN with small energy and momentum 
transfers attracts attention of experimenters as perspective tool for condensed matter studies 



[1J| . Recently observed (see, e.g., [14, [15|]) small cooling and heating of UCN in vessels also 
belong to quasielastic processes. 

A step from frozen to moving nuclei target is very dramatic since it requires a transition from 
one body to many-body function ^(r, R^). MST does not present any universal solution, since 
general equations of MST are only a reformulation of the problem in the way where multiple 
scattering is clearly exhibited (by using iterated (11) in (|To|)). Practical content of MST is, in 
fact, a set of approximations applicable for different situations. They were analyzed, e.g., in 
monograph of Goldberger and Watson || . 

Attempts were made to consider inelastic processes for UCN by using one of approaches de- 
veloped in MST (see, e.g., fl6|]). However, approximations used so far for MST cannot be applied 
for inelastic scattering of UCN. Indeed, the main assumptions, that different approximations of 
MST were based on, may be formulated as: 

(a) energy of projectile is much larger then characteristic energy of target particles ("weak 
coupling approximation" ) ; 

(b) mean free path of projectile in target is much larger then its wave length; 

(c) mean free path of projectile is much larger then length of effective correlation between 
target particles. 

The first of these assumptions allows to use Born or "impulse" approximation, where each 
target particle may be considered free when colliding with projectile. The second and the third 
condition allow to treat multiple scattering as sequential collisions and represent result as sum on 
number of collisions executed. Due to (b) the energy of projectile between successive collisions 
is well determined quantity, and due to (c) the target at each collision may be considered being 
in the ground state. 

It is easy to see that all three assumptions are not valid for UCN. 

(a) Energy of ultracold neutron is ~ 10~ 7 eV and corresponds to temperature ~ 1CU 3 K 
which is much smaller then the target temperature even at liquid helium. 

(b) The usual definition of mean free path fails for UCN (elastic cross section for mirror like 
potential is equal to surface area S). Thus, one may use for this quantity intrusion length into 
target, which is of the order of wave length ~ 10 nm. It means that neutron energy between 
collisions is, in fact, uncertain. 

(c) The main effect of UCN multiple scattering is appearance of a potential barrier, that is 
just the product of particle-particle correlation at distances compared with neutron wave length. 

The goal of this work is to find solution of MST equations valid for the whole domain of 
slow neutrons (from thermal to ultracold), starting from realistic neutron - nucleus interaction. 
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Equations ([|) and (||), (0) will follow from this theory as limiting cases. By solution we mean 
reduction of general equations dT0|)-(12) to those which allow reasonably simple numerical so- 
lution for elastic and inelastic scattering for all practically interesting cases. One numerical 
solution for inelastic UCN scattering is presented as an example in the last part of the paper. 



2 Formulation of the Problem. Plan of Solution 

A proper theory for UCN scattering should be based on the following postulates: (i) No Born 
approximation; (ii) No use of Fermi potential; (iii) Target matter is a dynamical system. So we 
should start from N+l body Schrodinger equation 

(j^+Ht + vj \^,i)=E Ki \V Ki ), V = Y^V V . (13) 

Here V u describes interaction of neutron with z^-th nucleus, p = —id/dr is the operator of 
neutron momentum, Ek,i = £k + £j is the total energy as the sum of neutron energy ek = k 2 /2m 
in the state |k) and the target initial energy e% in the state \i) that is the eigenstate of the target 
Hamiltonian Hf. Here and onward we keep fi = 1 till final physical results. 
Equation (^) can be written in integral form 

\^,i) = \K,i)+D- 1 J2Vu\^ k ,i), (14) 

V 

where l^w) = 1^)1*) arm D^ 1 is "Green function" with 

D = — + e i -^--H t + ir l , (15) 
2m 2m 

where positive quantity r] — > provides outgoing neutron wave asymptotic. 

Note, that the problem can be easily reduced to MST equations (|ic|)-(|l^). First, equation 
(14) can be written in the form 

|*k,i> - D^Vul^i) = |<*> + D- 1 J2 V^l*k,i>. (16) 

Then, let define a state vector and operator t u by the relations 

\^») = \^)-D- 1 V v \^ l ), (17) 

iu\*u) = V v \Vk,i). (18) 

Now one can see that equations (|l4] ) and (^) with the help of (|i~7|) and (jl^) turns into ([n]) and 
(11), respectively. Finally, we have to show that i v obeys (|j~2|). For this purpose we use identity 

ft|*k,i) = ft(|tfk,i) " D^V^i)) + VvD- 1 ^*^), (19) 
that, with the help of jl7|), (H), transforms into 

= v v \m v ) + v u d-%\^ u ). (20) 

Thus we have demonstrated that MST equations are nothing more than reformulation of the 
general scattering problem (|l^) or (|i~4|). 

It is, of course, impossible to solve the many-body equation ( |i~3| ) or MST equations (|i~0|)-(|l2|) 
and to found the state vectors |^k,i) or \^u) without any approximations. In our problem there 
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are two main small parameters: short-range of neutron - target nuclei interaction (as compared 
with interatomic distance and wave length) and small neutron energy (as compared with depth 
of interaction potential). 

The first condition allows to consider only s-wave part of the wave function of neutron - 
nucleus center-of-mass motion, when their interaction is evaluated. And the second condition 
allows in this evaluation to neglect energy of relative neutron - nucleus motion inside the 
interaction potential area. So, the s-wave function and its derivative, taken at the potential 
boundary, are independent on neutron energy and are just numerical parameters. 

No specific model for neutron - nucleus interaction potential will be needed. Its specific 
features described above (short range and large depth) allows to use scattering length approxi- 
mation. 

The small parameters allow to simplify our problem. Potential V v is essential only in small 
vicinity of R^. It differs from zero only when the absolute value of deviation x = r — R u does 
not exceed the potential radius r§ v . So, using completeness of the neutron states J2 r l r )( r l = 1 
one has 

vtl*k,i) = Ex|R- + x)(R, + x|v;|* kji ) = 

(21) 

= Ex |R* + x)K(x)(R^ + x|* k)i ). 

Here and onward sums of continue variables mean integrals with the following supposition for 
position- and momentum-energy variables 

XW(^ xWt 



Using (21) one can transform fll4| ) to 

l*k,i> = |* k)i ) + D- 1 E E l R - + x >K(x)(R, + x|* kji ). (23) 

V X 

The scalar product of (16) with (Rj, + x| after some rearrangement can be presented as 
(Rj, + x|*k,j) - 

- E( R ^ + XI^IR^ + X / )K(x')(Rz, + X'|* M ) = 

= <R, + x|*U+ (24) 

+ E E( R - + xl^ 1 ! 1 ^' + x , )K'(x / )(R.^' + x'|* M ). 

v'^v x' 

Equations ( p3| ) and (p^), as can be recognized by the structures of their right-hand sides, 



correspond to MST equations (10) and (|llj). To make them fully determined it remains to find 
the only key element, namely (R u + x^^i), i.e. the exact many-body wave function, but only 
in the area of short range potential for each nucleus. 

So, the many-body problem is reduced to only two-body problem (or, more precisely, three- 
body, since the nucleus is not free), with the rest of the nuclei as spectators. The solution of the 
last problem is determined by two parameters: by the values of s-wave amplitude x^( r Ov)/Vov 
and its derivative at the boundary of the potential. Since the logarithmic derivative is connected 
to the scattering length a u (the known physical quantity), the amplitude Xv{ r 0v) remains the 
only free parameter. 

Therefore one may hope to express the quantity (R^ + x|\Pk,i) by Xvfrou) and then to use 
(24) as a set of linear equations for parameters Xv( r 0v)- We will call (R u + x^^i) for x = ro 



V 



as " neutron function at the surface of u-th. nucleus" . 
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Remark. The small quantity r§ v will be neglected whenever possible, except in cases where 
Tq v stays near the scattering length (which may be of the same order of magnitude) and in 
terms with singularity l/rQ U till their compensation. Such singularity occurs in both terms 
in the left-hand side of (^4|) and requirement of their compensation will give us an additional 
control of calculations. 



3 Neutron function at the surface of z/-th nucleus 

Let us transform the basic equation (|l3|) to new variables 

r, R — ► x = r — R^, R, (25) 
where R = {R^}. The Hamiltonian in the new variables takes the form 



2m,, 



+ V u (x), 



(26) 



where p x = —id/dx, and P u = —id/dYl u is the momentum operator of u-ih nucleus. We 
assume that the quantity x is small, so that neutron interaction with other nuclei (with v' ^ v) 
is absent. 

We look for the solution of the Schrodinger equation with the Hamiltonian (|26| ) with the 
energy E^ i = k 2 /2m + £j. In Born approximation in the limit x — ► solution has the form 



*k,i(x,R) = ^(x)e 4kR ^(R) 



(27) 



where is the center-of-mass wave function, and $i(R) is the initial target state vector. In 
general case the neutron near the z^-th nucleus may have its energy different from the initial one 
due to previous collisions. So, it is natural to look for the solution in the form 



(28) 



where ^(R) is the eigenfunction of the target Hamiltonian Ht with the energy £j and g is a 
vector parameter which represents neutron momentum. 

If we substitute ( |28|) into Schrodinger equation with the Hamiltonian ( |26|) and take into 
account that operator P v acts on <&j(R) as well as on the exponent, then we obtain the following 
equation for (p(x) 



£r . g 2 . gGj„\ 
H n + TTT" + — rr- I V? 



(g + G ju )p x 



2M„ 



M,, 



V = (E^i - £j)tp, 



where 



G 



3 V 



(29) 
(30) 



Equation (E{J) after some formal transformation can be displayed as 



2m v 



in. 



+ V v {x) <p = E(R)<p, 



(31) 



where 



E(R) 



E u 



9 l 1 
— H 

2m 2m h 



g 



my 
Ml 



(g + G. 



(32) 
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Formally, ( |31| ) is an equation only for the function tp, and its exact solution, as can be easily 
proved, may be presented in the form 

p(x) = exp + G,v)xJ *(q,x), (33) 

where ^(q, x) is the scattering wave- function in the center-of-mass system determined by the 
equation 

tf n *(q,x) = ^-*(q,x), (34) 
and the scattering energy is defined from 



' + - — g - -rr(s + G 



2m u ' 2m 2m 



V 



(35) 



So, initial problem of neutron scattering on the bound nucleus seems to be reduced to the 
problem of scattering on free nucleus. In fact, it cannot be done precisely. Indeed, though 



(33) is formally an exact solution of the equation fl29|), its parameter q - impact momentum - 
through vector Gj u (|30| ) depends on all coordinates of the target R. But such a dependance 
was not assumed in (P8|) . 

However, we are looking for the solution valid for small x. In the limit of small x only the 
s-wave part Xv( x )/ X of the scattering function ^(q, x) is of importance and dependance of Xv 
on q can be neglected. Therefore, expression ([^) for small x is independent on R and gives the 
real solution in the form 

* k ,(x, R),^ ~ xMe**"*^). (36) 

So far the target state j and intermediate neutron momentum g were not specified. It is 
evident, that any linear combination of functions (36) is allowed (provided the right-hand side 
of ([?5|) is non negative), so we finally obtain 

\P kii (r = R„ + x,R)|^ ~ ^ e ^(R|^), (37) 

where 

<RW = *„(R) = £a(g, iK^^-CR). (38) 
s,j 

Note, that C u (g,j) = dg^Sji and Q u — > <&j in Born approximation. 

Expression ( |37| ) defines a local structure of the basic function near the point r = R„. 
It naturally contains a "background target function" 3V (R) = (R|^), separately defined for 
each nucleus. This function differs from the initial state of the target $i(R) = (R|i) due to 
perturbation by the neutron wave. The functions (R|^) should be consistently determined in 
parallel to the amplitudes of the neutron wave Xv( r 0u)- 

Let us now calculate the integral over x in (|2l|) . First, from ( |37"|) one finds 



(R^ + x|* k)i ) ~ ^M e ^|z,). (39) 

x 

Second, we note that the state vector \R V + x) due to smooth dependance on x may be factored 
out from the integral and taken at x = 0. Third, from the equation ([P|) for Xv inside the 
potential area it follows 

W ^ ( \ 1 d2 Xv(x) 

V u (x)xv(x) ~ 40 

2m v ax z 
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Thus, the integration can be performed 



(41) 



where we used the relation between scattering length a and the logarithmic derivative at the 
boundary of the potential [x'/x]r=r = — V( a ~~ r o)- Combining ( |39| ) and (fil|) we get 



V)|R I/ + x)V r I ,(x)(R v + x|* k)i ) = |R,)— e ikR »^, 

X 

where we have introduced the amplitude 



(42) 



(43) 



4 Equations for neutron amplitudes 

In the previous section integrals over x in (p4[) are expressed in terms of the amplitudes 4> v and 



state vectors \v). Thus, with the use of ( p9| ) the equation (24) takes the form 



x 1 ' 



(Pv + xlD-^R,,)— e* R » 



m 



»'W)(f>u>- 



(44) 



The small quantity x is left only in the terms with singularity 1/x and neglected elsewhere. 

Note, that x in fl44| ) is a free parameter, and for x = r 0u , when Xv{ r 0v) and <p u are linked 
by (f|3|), relation (||) may be regarded as a set of linear equations for <p v , but with operators in 
target states as coefficients could be hardly used for practical calculations. Our next steps are 
directed to simplification of (44). 

It is useful to transform the terms with as follows 



where 



D q = (q|D|q) 



q 

k 2 -q 2 



+ Si- H t + IT). 



(45) 



(46) 



"' ^' 2m 

Then from (44) it is evident that to separate the amplitudes one should multiply both sides by 



e jkR " and take their scalar product with the eigenvector (j\ of Ht- Thus, we have 







--) 


m 







where 



2tt 



C,W,x) = — V e iqx (jK i(k ~ q)R ^ ~V (k - q)R ^|z/>- 
m _ 



(47) 



(48) 



Expression ( |48| ) can be transformed into the other form with the use of the following pre- 
sentation of the operator D^ 1 



D- 1 = - I e i(^~^+e t -H t +i v )t dt 



(49) 
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Then we get the alternative expression for (|4 

2vr 



0(^',x) = ^x 

(50) 



oo 

X ^ 



^- e iqx j Xj (yy>^ k - q, t) e i{^-^i-e S +iv)t dt 



" 

Here the correlation function is introduced 

Xj (vv',K,t) = ( <7 -| e -i«R-(*) e i«R-„'(o)| z/ ') ) ( 51 ) 

which is a generalization of (0). 

To make (47) fully determined, we only need to find an explicit expression for diagonal term 
x) in the limit x — > r§ v . It evidently has a singularity l/rQ U which should compensate 
analogous term on the left-hand side of (flT]). That can be easily seen from simple scaling 
analysis. Let introduce new variables and parameters 

q' = qx, k' = kr, e\ — e'j = (ej — Ej)x 2 , t' = t/x 2 . (52) 

Then small parameter x remains in ( |50| ) only as the common factor 1/x and inside Xj i n the 
argument of TL u (t'x 2 ). 

When u = u' the operator in matrix element (51) is close to unity with very small deviation 
~ uk, where u is a nuclear shift from equilibrium. Thus, it seems to be a good approximation 
to integrate over t in (^0|) only the exponential factor outside the matrix element. Then the 
integration over q gives 

(j(w,x) ~ (j\v), (53) 

x 

where 

k 2 = k 2 + 2m(ei - Ej). (54) 

Note, that kj is the absolute value of momentum of inelastically scattered neutron when target 
remains in the eigenstate \j). 

However, this result is evidently wrong. Since it does not cancel the singularity 1/tq v in the 
left-hand side of (47) due to the additional factor m u /m. So, the matrix element Xj should be 
more accurately estimated. 

The limit x — » means that only small time interval is of importance in R^(t) = H u (t'x 2 ), 
and we can use the expansion 

K u (t) ~ R,(0) + ^. (55) 



Each additional power in t will result in an additional factor x 2 . After substitution of ( |55| ) into 
(|5l"|), the matrix element is calculated without any father approximations. Using the identity 
exp(A) exp(B) = exp(A + B + [A, B]/2), which holds when A and B commute with their 
commutator [A, B], we obtain 



— ► 01 exp - 



(k-q)P„ (k-q) 2 



\ , . (56) 
t\\u). 



Then integration over time and elementary transformation give 
(j(w,x)\x-*o = 

= iqX , \ ( 57 ) 

q 



n m 2^\J I _ + 2m v (ei - £j) + if] 1 
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where 

k„ = k-^(P„ + k), q„ = q-^(P„ + k). (58) 



The physical meaning of (|58|) is evident: \s. u and are the momenta in the center-of-mass 
system. This transition to the center-of-mass system provides additional factor m u /m. 
Finally, integrating over q and taking into account only s-wave part in x we get 

Tfl f 1 *• \ 

(j(w,x)\x->o = ( -UW) +i{j\ K juW) J , (59) 



where the operator Kj u , defined from 

+ 2m v {ei - Ej), (60) 



K 2 



k-^(p„ +k ) 



has the meaning of absolute value of the impact momentum in the center-of-mass system for 
the neutron and v-th target nucleus when the target is in the state \j). For UCN this quantity 
is defined mostly by the average absolute value of nuclear momentum. 

The limiting value (|59| ) allows to get explicit form of the equation ( p7[ ) for the amplitudes 
4> v . With the help of (|59|) we verify that for x = Tq v singular terms 1/vqv are canceled out and 
the next terms result in 

fl-OM^f + i—{j\Kjv\v)<pv = Sij + ^2 Cj{vv')(t) v <, (61) 

Pv m v ,^ u 

where the renormalized scattering length ("on bound nucleus") is introduced 

vn 

V = a v . (62) 

m v 

The matrix Q(i>v') = Q(uu',0) with the use of ( |50|) can be expressed in terms of Fourier 
transforms of correlation functions (|5l|) 

= 4vrV XiW±-w) , (63 ) 
w k j ~ 1 2muj + ir l 

Now we note that any state vector \v) can be presented as series in target eigenfunctions 

w)=J2\j)m- (64) 



Therefore the relations (31) are really the linear equations for the amplitudes 

€ = UWHu. (65) 

Indeed, we get 

l^ + ^^( i |^ v | i ')4' = ^+ J2 ^A""'H- (66) 

Pv Tfl .. . . . , 

]' 3 > v r v 

The coefficients of these equations 



q 

= l7T y X. 7 y(^,k-q,h;) 

J k'j — q 2 — 2muj + in 



(67) 



10 



with 

+00 

Xjj>{vv',K,u)= J Xjj>(w',K,t)e tujt dt, 

-00 (DO) 

are completely determined by the properties of target matter, i.e., by the matrix elements of 
the operator of v-th and is'-th. nuclei position correlation between the target eigenfunctions. 

The equations ( |66| ) are obtained from (p4|) or (11). However, they deal not with N+l body 
state vectors l^k.i) or \^f u ) but with numerical parameters <^ v . It is easy to see that due to 
(42) and ( |23| ) the total wave function l^k.i) can be expressed in terms of the same parameters, 
namely, 

l*k,i) = +D- 1 J2\K»)—j* R "\j)€- (69) 

It means that scattering probability and scattering cross section are also determined by (fy. 
Connection between them is analysed in the next section. 



5 Scattering problem and neutron amplitudes 

Scattering process with a fixed final state (j\ of the target is usually called a transition into j-th 
reaction channel. The wave function in the j-th. channel can naturally be defined by 

tftf(r) = (r,i|* k)i ) (70) 

Prom ( |69| ) it follows 

v^-(r) = 5 i3 ^ + £ ^(rlD-^y^lf). (71) 
Then using fl45|), (|4q) one has 

v^-(r) = %e to + 4tt ^ E k 2- r ~*+\/ kA ^> ( 72 ) 

j',v 1 3 " ' 



where kj was introduced by (|54j). After integration over q we finally obtain 

ik \r— R, I 

vE-(r) = <%e* kr 1 3 ^\ ^\f). (73) 

In asymptotic, when r — > 00, we get from ([/"3|) for the scattered wave the usual structure 

fijfakj)^, (74) 
where the scattering amplitudes in reaction channels are given by 

/«(k,kO = -E^'01 ei(k_k ' )A "b")- (75) 

The final momentum k' = k, is defined by kj = kj(r/r). 

It is instructive to note that one obtains in Born approximation with pseudopotential (|l]) 
the following expression for the scattering amplitude: 

^■(k,k') = -E^0> i(k ^ k,) ^K)- (76) 
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It is analogous to (|F5|) but instead of <ftl it has the scattering amplitude on an bound nucleus, 
and the matrix element is taken between unperturbed states (j\ and 

Remark. A typical target model used for ultracold neutron reflection is semi-infinite sub- 
stance. The asymptotic procedure used above to extract the scattering amplitude ( fT5| ) from 
the wave function, strictly speaking, is incorrect for an infinite target. So, we shall assume our 
target finite with some plane surface area S, which is large enough to use in parallel planes con- 
tinuum spectrum, orto-normalized with 5(ky — ky) instead of Sy^, but allows when necessary 
to make replacement 

"(2vr) 2 5(kj| - k||)] 2 — ► S(2iT) 2 5(k\ l - ky). (77) 

In the case when transmission is also of importance we need to consider the target of finite 
width to allow asymptotic in both directions. 

The cross section for the reaction i — > j can be obtained from the scattering amplitude by 



dai^^lfijiK^fdn', (78) 



but the final state of the target is usually not fixed and inelastic processes are measured by the 
energy transfer. This process is described by double differential cross section 

T E l^'( k > k ')! 2 ^ " ej + ek - £k0- (79) 



dU'de' k 

3 



It is useful to consider the quantity 



W (k,k') = £ 



m 



27r<5(£i - Ej + e k - e k /), (80) 



then w(k, k')dk'/(2vr) 3 is the scattering probability per unit time from the fixed state |k) to the 
states |k') into the momentum space element dk! . The scattering probability (|BC|) divided by 
the incident neutron flux k/m gives the cross section to find the final neutron with momentum 
k', and one is free to choose parameters in k' to be fixed in the final state and then integrate 
over the rest of these parameters. The relation between w(k, k') and the cross section is evident 
from the equality 

So, the cross section of inelastic scattering with energy loss uj is given by 

da m C . ,s r., s c?k' 



Replacing delta- function in (SO) by the integral over t and using ([75|) one can sum in <^ 
over the final states j as 



V'e^fe ~ £ j" )* (j | e -«(k-k' )R„ | j/> } (jii | e j(k-k')R,/ | j' } 
j" 

= Xjj>(w',k- k',t), 



(83) 



where Xjj'( uu 'i K it) ( as weu as its Fourier transform) was introduced by (|68|) . Then ( |80| ) takes 
the form 

w(k, k') = (—) Yl 4*<tt>Xjj> k-k'^ + Ei-Ej). (84) 
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The corresponding equation for double differential cross section is 
d 2 a k' 



— (^(/iajj'^u', k-k',u; + Si -Ej). (85) 



jjvv 

We have now the scattering amplitude ([75]) and scattering probability (|3|) expressed by the 
neutron amplitudes 4>l, which are determined by the equation (|66|). 

6 The case of thermal and cold neutrons 

Here we address to the neutrons with momenta in the range 



y/Aima < k < l/a, (86) 

for which the rescattering processes are not important and may be considered as small pertur- 
bation. 

In the first approximation the last term in (|6"6|) may be neglected and we obtain 



cpi = 5ij ^ — ^ 5ij(3 u (l - i{i\K iu \i)a u ) , (87) 

1 + i{i\K iv \r)a v v 7 

which is similar to the expression (||) for the amplitude of neutron scattering on the isolated 
nucleus. The only difference is that instead of impact momentum in (J2j) (natural parameter for 
two-body problem), in ( |87| ) enters the average value of this parameter over nuclear ensemble in 
the target. After substitution of (^7|) into ( j35[) we obtain the known expression (||) for thermal 
neutron scattering cross section. 



7 Renormalized amplitudes for condensed target 

Up to now no assumptions were made on the target matter. For a condensed target the results 
can be presented in more visual form. 
Let suppose 

H u = Pu + m p , (88) 

where for solid state target p v is the equilibrium position of the target nucleus, and u u is the 
shift from the equilibrium. For liquids p v and may be understood as an average position 
and a fluctuation. Note, that index v may be now replaced by p. When useful we shall make 
such replacement without notice. 

The matrix Xjj l { vv \ K ^) ls transformed now into 

Xj f(vu',K,t) = e - iK ^-P^(j\e- iK ^e iK ^%'). (89) 

For the model with fixed (unmovable) nuclei the matrix element in ( p9| ) equals to 5jji. In order 
to separate the effect of nuclei motion let us write 

Xjj '(^, t) = e~ iK ^~P^ (5 jf + Xjj'H, t)) , (90) 

where 

X jr (ui/,K,t) = (j| e -^M*) e ««v(0) _ il/). (9i) 

To simplify the equation ( |66|) let us introduce instead of Cjj'( w ') (|67| ) a new matrix Gjy(yv') 
from 

^,(yv') = -e~ i ^-P^G jr (yv'). (92) 
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Then introducing the new amplitude 



^) = i-4e ikp -', (93) 

Pu 

we obtain the equation 

V>» = dij^ - £ G3i'{w')4>j'{v')l3 v >, (94) 

j'u' 

where the diagonal in v and v' term is of the form: Gjj'{vu) = i(m v / rn)(j\Kj v \j') . Note, that 
according to ( |59| ) this diagonal term equals to — Cjj'( ul/ i x ) m the limit x — ► but without its 
real part in contrast with the non-diagonal terms. In the long wave length limit the sum in ( p4T ) 
may be replaced by the integral over v 1 . In this case it is not necessary to give special attention 
to the point v' = u, since the real part of the term with u' = v though singular but integrable 
and does not contribute into the integral. 

For the matrix Gjji(vv') we have from (|67[), ( |90[ ) and ( p2| ) 

Gjj'(yv') = bjj'Gjiyv') + Gjji(vv) (95) 

with 

Gj{vv') = -4vr£ = — — — — p (96) 

JJ V ; ^ k] -q 2 - 2muj + ir) v y 

where Fourier transform of Xjj'iyv' , K,t) ( |9~I| ) is introduced. 

The scattering amplitudes in reaction channels (ff^) in the new notations are given by 

fij^k') =-^e- ik '^0> i{k - k ' ) "1i')V'i'W^- (98) 

The scattering probability can be found with (|98| ) and (|8Q| ) or directly from (R4|) with (|9lf ) and 
©• 

Remark. In Appendix we obtain an alternative form of (^i|) and (|98|), which may be useful 
for some approximations. 



8 Unitarity condition 

The flux of neutrons incident on the target should be equal to the total flux in all scattering 
and reaction channels. As will be proved in this section, this unitarity condition is satisfied by 
the general theory suggested. 

It is convenient to rewrite general equation ( |94"|) for neutron amplitudes in the form 

^•e ik ^ = ^ (u) + G n' K) V)p\> , (99) 

j'u' 

and after complex conjugation 

%e-*^ = ^jiy) + E Gyiyv'Wyiy'Wu'. (100) 

j'u' 
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Multiplying (|99| ) by ip*(v)f3* and ( |100| ) by ipj(is)f3 u , then summing over j, u and subtracting 
one equation from the other we get 

2ilm/ ii (k,k) = 2(-2iImA,)^MV'iH + 

+ £ r j iy)$ j 'V)PtP V ' ( g jA™') - G *rM' v )) > 1101 ' 

where we take into account equation ( pg|) fo r scattering amplitude. 

The kernels G and G*, as seen from (|67|) and (^), differ only in path directions around the 
poles on complex q plane, and their difference can be presented in the form 

G jf (uu') - G*,>V) = -toLyeWu-Ps) x 

" " 171 q (102) 



x ^| g -i(k-q)u„ ^jj-l _ (£>-!) + j e Kk-q)V|y). 

Operator Z)" 1 has the form (a + irj)~ l and due to the equality (a + iry) -1 = — m5{a) + P(l/a) 
only delta- function part remains in combination D^ 1 — (D~ q 1 ) + . Thus ( |102j ) reduces to 

- G*,-(v'v) = dri'Ykjn^'tPv-P^ X 

■> j" (103) 

where n" = kj«/A;,«. 

Substituting this result in ( |101| ) and taking into account ( |98| ) we get total cross section 
as sum of capture cross section a c and total scattering cross section a s : 

a t = —ImfeQa, k) = a c + a s , (104) 

where 

- C = y(-^^M^\ 2 , (105) 

<t s = E|/ l^(k,k')| 2 dn'. (106) 

Note that capture cross section is determined by imaginary parts of scattering lengths a v and 
f3 u = ma u /m u (—Airlvaa/k is the capture cross section for free nucleus). 



9 Interim summary 

Let us sum up our results. To satisfy MST equations one has to find N state vectors ^ u and 
N operators t u . We have reduced the problem to the set of linear equations d66| ) for numerical 
parameters <ft v . The price for this reduction is additional index j which arises due to introduction 
of unknown state vector \u) for each target nucleus. However, if = {j\y)<j>v are found then 
one can deduce ^ v and t v . 

Indeed, equations (|B"I|) and ( |6"6| ) were directly obtained from (H), where the left-hand side 
has the meaning of (R^l^). Following back from (]6l| ) to (pi) one can deduce for ^ v 

= |R,)e 4k ^^- (l + ia v K v ) \v)(j> v , (107) 
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where 

Ku = Y^\j){j\K jv . (108) 

5 

On the other hand, if to compare ( |69[ ) with (|l0|), one can obtain combination 

t^|^> = _ |R„) e * R »<^. (109) 
m 

From (|10D and ( ^09|) it follows for t„ 

i = — |R,) e ik ^ fl + ia^) -1 ^e- ikft "(R v |. (110) 

This expression for t u is rather complicated due to non-commuting in target space operators 
K u and R,. 

Remark. In (|l~7| ) seems to be a state vector in N+l dimensional space (r, R). However, 
for short-range t u the total wave function j is really determined by *$> u at r = R^. This part, 
in fact, is given by (|107|) . 

For fixed nuclei K u — > k v is the center-of-mass impact momentum and from (|110|) follows 



<r|i>') - J v 6(r - K u )5(r - r% (111) 
m 1 + ia u k u 

that, with the help of (0) and (|62|), coincides with Fermi pseudopotential (|l|). Just in this 
approximation (||) and (0) follow from ([H]) and (|TT|), respectively. 

Fortunately, in general case all physical quantities can be obtained directly from the ampli- 
tudes (f>l or and matrices i v are not needed. Thus we reduce MST equations ([To| ) and 
(11) to systems ( |66|) for 0^ or (94) for 

Let us emphasize, that the only approximations used so far are those for the interaction 
potential. It was assumed a short range and relatively deep, what is equivalent to scattering 
length approximation for the interaction. 

A relation of our equations to the traditional theory for slow neutrons was demonstrated 
in section |(| For thermal and cold neutrons the structure of equations (^) may be radically 
simplified due to the physically justifiable neglect of rescattering of the secondary neutron waves 
in the target. 

Our main goal here is UCN. A relation to the traditional theory for the (elastic) scattering 
of UCN will be in details considered below. Physically evident, that transition to that theory 
should occur by the neglect of the thermal motion of the target nuclei, which provides the 
radical simplification of the correlation function (|39"|) : 



Xjj ,(vv>,K,t)~d jf e- iK (P»-P«'\ (112) 

The amplitudes of the thermal motion are indeed small as compared to the UCN wave length 
and the approximation ( |112j ) for Xjj' seams to be justified. But if the thermal motion is totally 
neglected it is no possibility for inelastic scattering. So, for inelastic processes we need to include 
the thermal motion but may hope for proper simplification since a perturbation procedure is 
justified. 

The rest of the paper is devoted to a perturbational solution of the main equation (|94|) for 
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10 Expansion over the amplitudes of thermal vibra- 
tions 



10.1 Zero order (elastic) approximation 

In the long wave length limit (ku <C 1) it follows from (pi]) 

Xjj'{uu',K,t) < 1, (113) 
which in its turn leads to the inequality 

Gj, y {yv') < Gj(yv'). (114) 
These inequalities open a possibility for the application of a perturbation theory. 



Physically ( j!13| ) and (114) mean that in this limit the basic process of neutron - target 
interaction is elastic scattering and the probability of inelastic processes is small. 

If we neglect Xjj' then we obtain from (^) an equation for the amplitude if)^' \v) in zero 
order approximation 

^f\v) = b^Pv - ]T G^i/)^V)A/. (115) 



As seen from ( |115| ), the equations for different j are separated, and since inhomogeneous term 

4>f\v) = 5 lJ i>{v), (116) 



contains the factor <5jj, we have 



where is defined by the equation 

= e *P* _ £ G{yv')^{v')p v , (117) 



with the matrix 



G^)^G,^) = -A^ W - ¥T - = J —— V (118) 



Note, that singular real part of diagonal (u 1 = v) term in ( |118| ) as discussed after ( [94] ) should 
be extracted that results in G{vv) = ik. 

The equation ( 117 ) is similar to formula (|9|) customary used for UCN in the target model 
with fixed (i.e. in fact, infinitely heavy) nuclei and transforms to it by redefinition x I / k( l/ ) = 
ip{v){l + ik[3 u ) [12]. Its solution is basically simplified when the sum over v is replaced by the 
integral over p. Then after acting on (117) with operator A + k 2 it is reduced to Schrodinger 
equation with the potential 

2tt 

U = —n(3, (119) 
m 

where the density of the target n and scattering length (3 may depend on p. 

Formally, the equation ( |117| ) after replacement of the sum by integral J2u ~^ I n dp, becomes 
of a linear integral type. Such a replacement, if useful, will be performed below without special 
notice. On the other hand, to make presentation of general formulae in more compact and 
transparent form it is useful to consider all G{vv') as matrices, ^(v) and e* k ^" - as columns, 
^(v) and e -lk ^" - as rows and omit summation indices v and v 1 . In this notation (|117| ) looks 
as 

$ + = e ik ^. (120) 
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As seen from (118), the kernel G(k, vis 1 ) depends on the absolute value k, but a solution ^>(k, v) 
depends on the vector k, defined by inhomogeneous term on the right-hand side. 
The scattering amplitude ( |98| ) in zero order approximation is given by 



J^(k,k') = -M(k,k') 



:i2ii 



where special notation is introduced for /3-weighted Fourier-transform of the amplitude ^(k, v) 

^(k, q) = J] e~ iq ^V(k, u)/3u. (122) 



Scattering probability fl80| ) in zero order approximation is of the form 



(°)(k,k') 



Thus from (^) we obtain differential cross section of elastic scattering 



where Ik'l = Ikl. 



(123) 



(124) 



10.2 Approximations of the first and the second order. Inelastic 
scattering 



Zero order amplitude given by (121) corresponds to elastic scattering. Thus, the probability 



(80) for inelastic scattering starts from the second order term which is determined by the first 
order scattering amplitude 



^(k,k') = 2.£ ^/£)(k,k' 



5{ei - Ej + e k - e k ' 



Introducing an expansion in ku for neutron amplitudes 



^•M = «( J #) + ^ ) W + ... 1 
/«(k,k') = -^ (1) (k,k')- 

-i{k a -k^)Y^e-^'P^j\uZ\i)^)Pv 



one obtains from (98) 



(125) 



(126) 



(127) 



(i) 



where the first term is /3-weighted Fourier component (122) of the first order amplitudes ipj 
defined by the equation 

(128) 



Matrix G { 3 (yv r ) is the first order term which follows from (|9l| ) and ( p7[ ) 



jj 



Gf}{vv') = (V- - ik?) ((j\K\j')GAur/) - 

-G j {vp , ){j\ul,\j')), 



(129) 



where V£ = d/dp a v . 
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Fourier component ip^ 1 (k, k') can be found without an explicit solution of equation (|128|). 



3 

First, let us introduce the function 

$(k» = ^(-k», (130) 

determined by the equation 

^ + p^G = e' ik 'P, (131) 



where we consider tjj as row, and G{vv') = Gj(yv r ). Then, multiplying (128) by fUtf from the 
left and using ( 131[ ), we obtain in the left-hand side just the Fourier component that we are 
looking for. Therefore we have 

^ 1) (k,k') = -/3^ ) ^- (132) 



Right-hand side of this equation can be simplified with the help of ( |120| ), (131) and ( |129| ). 
Then, we finally obtain 



vf^k') = ^O'KH)v^(^WM) 



(133) 



V 

Inserting (|133| ) into ( |127D we get for the first order scattering amplitude 

4 1} (k,k) = - Y,M\u u \i)V u ^(k',^(k^)) . (134) 



Remark. An alternative derivation of expressions ( |121| ) and ( |134| ) is given in Appendix. 
The probability of inelastic scattering (in the second order) follows from (|125|) and ( 134 ) 



4 2) (k,k') = (^) 

x V- (Vi(k, i/)^(k, u))* VI, (^(k', z/)^(k, I/')) x 



(135) 

+0O 

x f e^-ej+^-e^t (i\u°\j)(j\u T u ,\i) dt. 



— oo 



Summation over j can be explicitly performed according to 

^e^-^^Ki^OK'K) = <»K(t)<&(0)|i>. (136) 
i 

This diagonal matrix element may exhibit a spatial dependence only as a function of — p u , 
and therefore allows a Fourier transform 

Then, we finally obtain for the probability of inelastic scattering 
ui 2) (k,k') = 

= - w)5**(q)iT(q)fi<-(q,a;), ' ''^ 
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where 

B(q) = XXe-^-V, $(k>W(k,z/)) . (139) 
The cross section for neutron to lose energy ui can be calculated from (^) and ( |138| ) 
^ (2-) 2 



duj mk 

q,k' 



E 5 ( e k-ek'-w)^(q)i? r (q)^ CTT (q,a;). (140) 



To disclose physical meaning of ( |140| ) it is instructive to compare it with corresponding ex- 
pression which followed from (||) . To make the comparison one should transform the correlation 
function x(j/z/,k,u) to f2 CTT (q, uj) (|137|). Expanding in ftu by the use of (|88[) and 



(H) we get 

x(»/i/,k,w) = 2-K8{uj)e- iK ^-Pw) (1 - (( KU „) 2 )) + 

T ( 141 ) 

+ e -*K{p v -p v >) l tf R T / (i|^(t)^,(0)|i) e ^*dt. 



Now it is easy to see that cross section with energy loss (see ( plj) and (p^)), followed from 
@ and flUD, 

can be reduced to the form ( p. 40]) with 

B(q) = i(k - k') ^ bve-WveiQt-^P* (142) 
i/ 

instead of B(q) (|139| ). The difference (apart from factor N) is that functions ip(k, v) and V>(k', ^) 
in (139) are replaced by the plane waves e* k ^ and e _lk P", respectively. It is very natural since 
(|3|) is obtained in Born approximation. 

So, one may say that ( |140| ), similar to (0), takes into account interference of two scattered 
waves (that result in inelasticity), but in addition uses wave functions in both input and output 
channels modified by rescattering. 

The idea to modify (^) for UCN by replacing plane waves with solutions of equation (^) is 
very natural and was tried in several papers (see, e.g., fl7|| ). But it is evident that if to do it in 
( |142j ) the result will not coincide with ( 139|) . 



11 Choice of correlation function 

Inelastic processes with energy and momentum exchange between neutron and target are es- 
sentially determined by the dynamical properties of target matter, i.e. collective excitations 
that are suitable (for given conservation lows) to provide this exchange. Correlation function 
that enters into cross section just describes these dynamical properties. The physical meaning 
of correlation functions - description of space-time evolution of a fluctuation appeared at some 
moment in some position point. 

The field of correlation functions is covered in a number of books and review articles (see, 
e.g., [||, 18, p^fl). So, here we'll just mention a few details necessary for what follows. 



Our function (p.37|) is related to density fluctuations 

2 ^(/n(r' + r,t)n(r',0)dr') =E eJqM ^(q^)- ( 143 ) 



Fourier transform S(c[, uj) (often denoted as "dynamical structure factor") can be shown to be 
connected with ( |137| ) by 

ft 

stf,u)~— wnrfau), (144) 
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where the quantities q' and q are equal but for crystals may differ by a reciprocal lattice vector. 
For simple model of harmonic crystal one can easily obtain (for phonon occupation factors 



where T is temperature and s is velocity of sound. 

Sound branch of excitation is effective for large energy and momentum transfer (say from 
UCN to thermal), but very ineffective for small transfer (of the order of initial energy and 
momentum of UCN). The reason is in fact that neutron dispersion low e ~ vk is quite different 
from that of sound uj = sq because for UCN v/s ~ 10~ 3 and one cannot satisfy two requirements 
Ae ~ uj and Ak ~ q simultaneously. For small transfers we need excitations with small uj and 

q- 

Limiting value of correlation function for u,q —* is given by " hydrodynamic value" 

n^(q,u) ~ 2 a ® . , (146) 

v ' q 2 nMs 2 uJ 2 +D 2 q i v ; 

where D is coefficient of any diffusion- like process, e.g., self- or thermo-diffusion coefficient (in 
the last case a = cp / 1 cy — 1, in the first case a = 1). At normal temperature parameter 
a = cp/cy — 1 is of the scale of 10~ 2 for solids and of 10 _1 for liquids. 

Function (|146j ) for fixed q has a pick value for uj = and width ~ Dq 2 in contrast to (|145| ), 
where uj and q are strongly coupled (uj = sq). It is useful to introduce instead D a dimentionless 
parameter 

2mD 

«* — s— ■ (147> 

which appears if one considers dimentionless variables uj/e^ and q/k. One may expect that 
the optimal conditions for small energy e ~ fiuj and momentum k ~ q transfer would be when 
this parameter d is of the scale of unity. In reality at normal temperature d varies from ~ 10 3 
(metals with high thermoconductivity) to ~ 10~ 2 (self-diffusion in liquids). 



The total correlation function include the phonon part (145) as well as all types of diffusion- 



like parts ( 146 ). All these parts are linearly summed in cross section and their contributions 



may be calculated separately. 



12 Subbarier inelastic scattering 

To consider a specific inelastic scattering problem with general formula (140) one needs, first, to 



find zero order (elastic) neutron amplitudes for input and output channels ^(k, v) and V>(k', v) 
and, second, to choose a correlation function that is adequately describes collective excitations 
of target matter in energy-momentum domain of interest. 

We consider, as illustrative example, scattering on thick uniform plane target when neutron 
energies in input and output channels are both below potential barrier. Let z axis is perpendic- 
ular to the surface of the target located at z > 0. After replacing discreet variable v by uniform 



p, one can reduce integral equations ( |120D and ( |131|) to Schrodinger equation 



(k 2 + A)V>(r) = u(z)iP(r), (148) 
where potential u(z) = 47r/?n(z) is determined by the target density 

»(*) = {"' HI (149) 
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Equations for the cross section contain the values of an d V>(r) only inside the target, but 
solutions for them are determined by the inhomogeneous terms of integral equations ( |120| ) and 
(|131| ) which for solutions of Schrodinger equation ( |148| ) have the meaning of waves incident on 
the target. We expand neutron momentum k in input channel into components ky and k±e z 
along and normal to the target surface. We use in what follows that k\ < k 2 < uq = A-K(3n. 
The solution of (|148|) for incident neutron in the region z > is of the form 



V;(r)=te ik H r H- 3ez , t=-^-, ae = J«o - (150) 

k± + tee v 

Neutron momentum in output channel we also expand in sum of longitudinal and transverse 
components k' = kf, + k'j_e z . Thus for the elastic scattering amplitude (121) we get 



/(°)(k,k') = -^(k,k') = 27Tik ± 5( 2 \k\ l - k„)^±*E (151) 

11 11 k± + Z8S 



where k' ± = k± for transmission and k' ± = —k± for reflection. Note, that here we take into 
account diffraction forward scattering (k' = k) which results from finite transverse size of a 
target (see remark related to equation ([T~ 



Substituting (|151| ) into ( 124 ) and integrating over solid angle around the direction k' = 
k|| — k±e z , we obtain zero order cross section of neutron elastic reflection from the semi-infinite 
target 

4 0) - 4- < 152 > 

Here S is area of plain target surface, and k±/k = cos 8, where 9 is the angle of incidence. That 
is simply the whole target area seen from incident neutron direction. This result is natural for 
total reflection. 

To calculate inelastic scattering to neutron state with momentum k' we need the solution 
ip(k' ,r) = k', r). Since subbarier neutron in output channel is back scattered it is convenient 
to assume that k' = k| — k' ± e z , where k' ± > 0. Thus 

2k' 



{[,(v)=t'e-^- 3i '\ t'= /, 7 , ^ = Ju -k'l (153) 

k I + ise' v 



Vector B (|13S|) is given by 



B(q) = n/?(2vr) 2 ^ 2 )( k|| - kj, - q„ )tt' ^L _ ^ + ^ . (154) 

For simplicity we neglect imaginary parts of amplitude (3 and potential Uq related to radiative 
capture. Thus, when substituting ( |154|) into ( |14C| ) we will use 

ft =W^-2=\ ( 155 ) 

k | + se 2 imp 



and the same for t'*t' . Taking into account (77), one gets for inelastic cross section the following 
expression 

(2) 

(l(7 ' : S^-rr^^h [dk'6(k' 2 + 2muj-k 2 )x 
k nMs z 7r 4 M J 

1 1 (156) 



duj 

xk' 2 /dq<5^(k|| -k' - qn) A(*4,q). 
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Here 

( o uj 2 \ <? 2 + (ae + ae"i 2 
A^(A4,q)=vr5 U 2 -^ 4 I I 



s 2 ; g 2 + (ae + ae') 2 (15?) 



for phonon correlation function (|145|) and 



A!*M a) ~ f g » +(3e + ae/)2 - f\ (158) 

for hydrodynamic one (|146[ ), where T = ^\uj\/D. Dividing the inelastic cross section over the 
transverse target area Sk±/k, we get differential probability per one bounce dwi e /duo for neutron 
transition to the state with the energy e' = e — lo. 

Integration in (156) over qii removes two-dimensional delta-function. Then remaining delta- 
function allows to perform integration over M, . The result is 



\/k 2 — 2mw 

dwi e _ T k± 
duj ~ nMs' 2 tt 4 \lo\ 



+oo 



k dkj^ x 
x J dip J dq ± A(k' ± ,q), 



(159) 



where one should keep in mind relations 

q 2 = k 2 + fc(| 2 - 2fc|| jfefj cos ip, k\ 2 + kf = k 2 - 2mw. (160) 

For the phonon model one can easily proceed further using delta-function in A ph and the 
small value of q 2 = q 2 + q\ = uj 2 / s 2 -C k 2 . Integration gives 



duj ~ 7T Ms 2 U s V C7 



(161) 



where &o = y2mU/h is the momentum at the potential barrier C7, v± = vcosO is the normal 
component of the incident neutron velocity, and ey = e sin 2 6 is the energy related to the incident 
neutron motion along the surface plane. 

For the hydrodynamic model we can perform in ( |159j ) integration over q± by closing inte- 
gration path in complex q± plane and suming over pole residues. Then it is useful to write the 
result in the form similar to ( |161| ) 

dw ie d _ 2 aT k ^ f (p n J ft 

where dimensionless function /, depending on e, 6, e' and parameter d = 2mD/h, is given by 



\/k 2 — 2'muj 

f(e,6,e',d) = — / k' 2 dk' ± I dyL(k' ± ,<p). (163) 



nd vq 

u 

Here vq = hko/m is the boundary neutron velocity, and 

T(U > >/2 (1 - 4A 2 ) 3 / 4 



1 + - 



1\ /i 2 + A 2 A(A + 1) 



H ) 1 + 2/i + 2(/i 2 - A 2 ) 1 + 2A 
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(164) 



Figure 1: Function f(e,8,e',d) from ( |163|) for fixed initial neutron energy e = U/2 and angle 
incidence 9 = 7r/4 versus final neutron energy e' and parameter d. 



2A 2 = j 2^ = ( " + "^ . (165) 



/\2 



To demonstrate the magnitude of the function / ( |163| ) and its dependance on the final neutron 
energy e' and parameter d we have performed numerical calculation of double integral for typical 
values of initial neutron energy e = U/2 and angle of incidence 6 = it/ '4. The results are shown 
by solid lines in Fig.|l| for d -C 1 and in Fig.^ for d > 0.1. 

It is seen that the spectrum of inelastically scattered neutrons has a peak in the vicinity of 
the initial neutron energy. In this region parameter T = y/\u\/D = ko^/\e' — e\/(U d) is small 
with respect to ae + as' and the main contribution into the integral comes from small q^(k' ± , (p). 
Then one may simplify the function L (|164|) by taking limit ji 2 3> (A 2 , 1): 



1 (1 - 4A 2 ) 3 / 4 
L--^^- — ^-4 — , (166) 
v^r 3 1 + 2A ' K J 



which allows to evaluate in (|163[) all parameter dependence and obtain 



Here C = 0.47 is the value of a dimensionless integral 



f(e,e,e',d)^C^^ A l—^\. (167) 



Approximation (167) is valid for |e' — e| = \u\ <C Ud, and for large parameter d ^> 1 ( |167| ) 
may give a good estimate not only for the small to peak but for the whole subbarrier area (as 
seen from Fig.^). 

Dependence on to is governed mostly by the parameter d. For d <€i 1 the peak is more 
pronounced and, when d decreases, becomes more narrow and high but with fixed (independent 
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Figure 2: Function f(e,6,e',d) for fixed e = U/2 and 9 = 7r/4 versus final neutron energy e' and 
parameter d. Solid lines - result of exact numerical calculation ( |163| ), dash lines - approximation 



on d) area. Indeed, using approximation (167) one obtains 



w 



hyd 



dw, 



hyd 



doo 



-dw 



Ud 

iC aT kjj 1 

vr ~M? U J V ^ 
o 

8C aT , fe± 
— , r '> k±p \ -ft. 
vr Ms 1 ±y yU 



duj 



(168) 



Contribution to the peak from phonon model can be neglected since it has smooth behaviour 
and contains a suppression factor v/s ~ 10~ 3 . 

For d > 1 the small u> peak becomes less pronounced (see FigfJ) and the probability for 
neutron to remain under the barrier after inelastic scattering diminishes as 1/ y/d (since in ( |168| ) 
the upper limit of the integral is now U). 



13 Conclusion 

The general theory of neutron scattering is presented, valid for the whole domain of slow neu- 
trons from thermal to ultracold. For thermal and cold neutrons, when the multiple scattering 
in the target can be neglected, the cross section is reduced to that known for thermal neutrons, 
which is determined mostly by correlation function for the target matter (section |6|) . 

For UCN the rescattering is the dominant process, but the theory can be simplified by 
exploiting small parameter ku, i.e. the ratio of the amplitude of thermal vibrations (for solid 
targets) or relaxation lengths (for liquids) to neutron wave length. In zero order approximation 
in ku (that is equivalent to the scattering on a target with infinitely heavy unmovable nuclei) it 
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follows the known equation for elastic scattering of UCN . Dynamical processes in the target 
are taken into account in the next orders in ku and result in inelastic scattering. 



A detail analysis of inelastic scattering needs separate publication. Here in section 12 a 
specific example was considered: scattering with small energy transfer when scattered neutron 
remains below potential barrier. This quantitative example allows to make some conclusions. 

The value of cross section is very sensitive to correlation function used. Phonon model 
which gives main contribution for UCN excitation into thermal region is quite ineffective for 
small energy transfer when space-time correlation processes are determined mostly by relaxation 
(" hydrodynamic" ) processes. So, the first condition to obtain reasonable theoretical result for 
cross section with small energy transfer is the choice of an adequate correlation function. 

The second factor that needs a reasonable physical modeling is elastic potential suited for 
target matter in each specific experiment. Consistence, uniformity, possible existence of surface 
layers, presence of hydrogen and its distribution - all that may require a change of model 
potential and therefore the wave functions in input and output channels that enter the cross 
section. 

All effects connected with neutron spin are outside of the scope of this work. For target 
nuclei with non zero spin the scattering length depends on spin-spin orientation. This would 
require only replacement in all formulae of scattering length [3 by weighted average value. (The 
same prescription is valid for isotope non- uniform target). For large wave length of UCN the 
averaging is well justified. 

Neutron spin interaction with target electrons (" magnetic scattering" ) does not present any 
specific difficulty, but require inclusion of new "spin-spin" correlation function. 

Quite different feature have spin-flip processes. Physically interesting is calculation of depo- 



larization probability for stored polarized UCN. This effect (as well as neutron capture (105)) 
belongs to incoherent processes. They are not considered in this work. It is evident that in- 
coherent processes can be considered by simple perturbation theory with "elastic" functions as 
zero approximation. 
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Appendix: An alternative expression for scattering amplitude 

We start with transformation of equation (p4j). First, using the definition (p2]) for the matrix 



Gjji{vv') and the expression (67) for the matrix Qji^uv') one obtains for Gjji(vv' 



m, t-^i 



m V (ai) 



x (jle^O-qW £)-i e «(k-q)iv 



Then it is convenient to introduce the operators 

PjA v ) = (j\e-^+^^\f), 



(A2) 

P jj >(,v) = (j\e i ^- iV ^\f), 

where arrows on operators Pjji{v) denote direction for gradients to act on functions of p v . Thus, 
after evident transformation with the help of (|96| ) the matrix Gjj'(vv') takes the form 

Gjji(vv') = P jjii{y)Gjn{vv')P ji'j'iy')- (^43) 
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Equation (|94|) with the use of (A3) can be written as 

^•0) + E p jj»{v)G r ,{vu')p j // i /(z/)VyO / )A' 



(A4) 



(A5) 



Note the action of operators Pjj> (v) on exponents with k and k' 

e - ik 'P»P~ jf (v) = e- ik 'P»(j\e i( - k - k >»\j'}. 
The scattering amplitude (|9|), due to (A5), can be represented in the form 

filiKk') = -J2e~ tk 'P»P jf (»)iP r (v)f3 u . (Ad) 



TV 



Now let act on all terms of (A4) by the operator P~}-{y) and sum over j. The right-hand 
side, due to (A5), remains unchanged and we obtain a new form of general equation (p4]) where 
the matrix Gj(yv') is "open" from the left 



ji ji u i 



(A7) 



We may now multiply (A7) from the left by /^i/^ (^)> which is a solution of equation ( |131| ), i.e. 

Y,Mf\v)G 3 {w') = e" ik '^' - $V). {AS) 
Then summing over v we get with the help of (A8) 

Y,Mf\v)P^{y)^{y) + 

j'u 

V 

From (A6) and (A9) it follows for the scattering amplitude 

f ij (k,k') = -s ij j2^^4 0) ^)- 

V 

j'u 

where zero order term is explicitly extracted. 



(A10) 
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As the last step, we perform the action of the operators Pjj/(v) on $j (y) and ipj(v) and 
arrive at desired relation between scattering amplitude and exact solution ipj (v) of the general 
equation (|9~4|) 

/ ii (k,k') = -%E^ eikP ^f ) ( i/ )- 



J2(j\^ u "{-<pf\p» + ^)Mp»)- ( Aii) 

j'u 

-i>f\p v )Mp v -^u))M). 



An expansion of (All) in u u gives in zero order (with the use of ( |130| ) and ( |122|) ) 

4 0) (k,k') = -^V(-k , ,-k), (A12) 

what, due to time reversal invariance, equals to ( |1 2 1[) . The first order term in (All) coincides 
with (1133). 
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